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^ ■ 1 Introduction 

■ The formation of steady spatial structures in systems far from equilibrium has been studied in 

j^ . great detail over the past years, the classical examples being Rayleigh-Benard convection and 

Tayor vortex flow ||^]. In quasi-one-dimensional geometries they usually share a common feature: 

'^ ' the stable structures that arise after the decay of transients are strictly periodic in space and if 

they are weakly perturbed they relax diffusively back to the periodic state. This relaxation has 



been investigated experimentally in various systems [E3] and the results agree with theoretical 



C^ 



*^ ' results based on the phase diffusion equation |21, |13|, |15(|. It describes the slow dynamics of 



f— ^ ■ the local phase, the gradient of which is the local wave number. The structures are generally 

^^ , stable over a range of wave numbers which is limited (at least) by the Eckhaus instability which 

^J I is characterized by a vanishing of the diffusion coefficient. It leads through the destruction (or 

T-j- ■ creation) of one or more unit cells, i.e. a roll pair in convection or a vortex pair in Taylor vortex 

CK , flow, to a new stable periodic state with a different wave number. 

r— H I Recently it has been pointed out that a vanishing of the diffusion coefficient does not neces- 

^ • sarily invoke the Eckhaus instability Q, ^]. Instead, under certain conditions, the system can 

go to a state consisting of distinct domains in which the wave number has different values and 
C^ I which does not relax to a strictly periodic structure. The domains are separated by domain 

^H' walls in which the wave number changes rapidly. This situation has been treated successfully 

using a higher-order phase equation 0, p^j . 
k> , Experimentally, inhomogeneous structures have been found in a variety of systems. Most of 

Vh ' them involve at least one time-dependent structure, e.g. counterpropagating traveling waves (or 

spirals) |l^, (turbulent) twist vortices amidst regular Taylor vortices Q, localized traveling- wave 
pulses in binary-mixture convection [g] and steady Turing patterns amidst chemical traveling 
waves Ip^ . Recently, however, domain structures involving only steady convection rolls of two 
different sizes have been observed in Rayleigh-Benard convection in a very narrow channel p|. 
It has been speculated that these structures may be related to the phase-diffusion mechanism 
discussed above Q. 

Within the framework of the phase equation domain structures consisting of an array of 
domain walls are not stable due to the attractive interaction between the walls. This leads to a 
coarsening of the structures during which domain walls annihilate each other. This dynamics is 
closely related to that observed in spinodal decomposition of binary mixtures after quenches into 



the miscibility gap |1C]. In general, the coarsening will eventually lead to periodic structures 
with a constant wave number. Only if the boundary conditions conserve the total phase, i.e. 
require the total number of convection rolls, say, to be constant, will the final state consist of a 



(single) pair of domain walls. 

Here we investigate the possibility of stable domain structures in the absence of phase con- 
servation. To do so we study systems in which two periodic patterns differing only in their wave 
numbers are equally likely to arise. In certain cases the competition between the two wave num- 
bers can be described by an extended Ginzburg-Landau equation. We study this equation with 
spatially ramped control parameter in order to allow the total phase to change and find that 
even in this general case domain structures of varying sizes can be stable. This Ginzburg-Landau 
equation can also be viewed as a one-dimensional version of the Ginzburg-Landau equation for 
two-dimensional patterns in isotropic and anisotropic systems. The domain structures corre- 
spond then to 'zig-zag' patterns. In two dimensions phase conservation is less common than in 
one dimension due to the possibility of focus singularities in the pattern which often arise at the 
boundaries @ . Our results may therefore also shed some additional light on the behavior of the 
'zig-zag' patterns studied previously |17, ^. 

2 Numerical Simulation of Fourth- Order Ginzburg-Landau Equa- 
tion 

A situation in which the competition between two wave numbers can be investigated with relative 
ease is obtained if the neutral curve which marks the stability boundary of the basic state has two 
almost equal minima at (slightly) different wave numbers. For values of the control parameter 
T, slightly above these minima the pattern can be described by a Ginzburg-Landau equation for 
the amplitude A |18], 



drA = D2d\A + iD^d\A - d\A + T.A- \A\'^A, (1) 

which gives for instance the vertical fluid velocity in convection via 

v^{x, z, t) = ee^'?=^.4(X, T)f{z) + h.o.t. + c.c. (2) 

For simplicity we assume in the following that the neutral curve has reflection symmetry with 
respect to qc and set D^ = 0. In this case the neutral curve has two minima if D2 is negative. 
The control parameter S is proportional to the temperature difference across the fluid layer, say. 
In order to allow the number of rolls to change freely without any pinning by the boundaries 
we apply a subcritical spatial ramping to S. Thus A goes to zero before the boundary of the 
system and the wave number selected by the ramp |11, ^] corresponds to one of the minima of 



the neutral curve. For the numerical simulation a Crank-Nicholson scheme with dx ~ 0.034 is 
used. Such a small grid spacing is required to reduce the pinning by the grid below the small 
attractive force between domain walls. 

For slowly varying wave numbers Q(t, ^) = 9j<^, the Ginzburg-Landau equation (|l|) can be 
reduced to an equation for the phase of the amplitude A = Re^'^, 

d^cj) = {D + Ed^ct) + F{d^(t)f)dl(t) - Gdl>p. (3) 

Domain structures arise for negative values of D. For large values of S wave-number gradients 
across domain walls are sufficiently small to be described by (y), and within this framework the 
interaction between domain walls is purely attractive. Thus, in the presence of a subcritical 
ramp domain structures are expected to evolve to a strictly periodic pattern. This is confirmed 
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Figure 1: Evolution of the total phase without phase conservation. The corresponding final 
states are shown in fig.0. 



for T, = 100 and D2 = — 1, using initial conditions with three domains, i.e. a region of low wave 
number between two regions of high wave number. 

For S = 1, on the other hand, the behavior is quite different, as shown in figjl|. It gives the 
total phase / Q dX in the system (between the ramps) as a function of time for different initial 
conditions which are obtained by changing the width of the central (low wave number) domain. 
While for large T, the system evolves to the same periodic state independent of the initial total 
phase, it converges here to one of several final states with different total phase. These states are 
pictured in fig.0. They differ in the width of the central domain and, most strikingly, their local 
wave number exhibits oscillatory behavior in space. 

The possibility for spatial oscillations in the wave number can be seen in a linear stability 
analysis of the periodic state. It shows that solutions which approach a periodic state for 
X ^ ±00 can do so in an oscillatory manner if the fourth derivative is present. The purely 
attractive interaction between domain walls within the phase equation, and therefore also for 
large S, is due to the monotonic behavior of the wave number across the wall. The non- 
monotonic behavior found for smaller S therefore strongly suggests an oscillatory contribution 
to the interaction which would explain the stability of bound pairs of domain walls. 

To study the transition between the two regimes one could investigate the persistence of the 
domain structures when increasing S. Since the scale for T, is set by D2, which determines the 
depth of the wells in the neutral curve, increasing T, is equivalent to making D2 less negative. 
In fact, one of the 2 parameters could be scaled away (L>2 -^ D2/\T,\,'E, -^ ±1). We therefore 
investigate the persistence of bound pairs by changing D2 rather than S since this also sheds 
some light on the stability of the two-dimensional patterns discussed below. The result is shown 
in fig.H. It gives the total phase within the unramped region of various states as a function of 
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Figure 2: Coexisting stable final states obtained from the time evolution shown in fig.^. The 
thick line gives the local wave number in the bulk (800 < i < 2400, dx = 0.034). Note that the 
amplitude goes to zero toward the boundaries due to the subcritical ramping. 



D2. For these simulations the system has been chosen much larger in order to avoid interactions 
between the domain walls and the ramps. The periodic state is indicated by a dashed line. The 
different symbols denote where the bound pairs with 1 to 5 minima in the local wave number, 
respectively, disappear. While the states with 1 to 3 minima jump to the periodic solution 
above that value of D2, the 4-state jumps to the 3-state and the 5-state merges with the 4-state. 
The 5-state is relatively hard to track due to the rapid decay of the oscillations away from the 
domain walls. For D2 < —0.8 all six states coexist stably. In this regime arrays of domain walls 
should be possible in which the widths of successive domains alternate chaotically Q. Clearly, 
none of the bound pairs investigated persists all the way to D2 = 0. This corresponds to the 
previous result that for large S the phase equation becomes valid. Strikingly, the value of D2 at 
which the solutions cease to exist is not a monotonic function of the number of minima in the 
wave number. Rather, the solution branch corresponding to the state with 3 minima is the last 
non-periodic state to disappear. If one were to start with a chaotic array of domain walls and 
to increase D2, fig.(^) suggests that domains of various widths would successively be eliminated 
beginning with the longest and the shortest ones until only domains with 3 minima in the wave 
number remain. Note that these simulations only address the interaction between two domain 



walls. In a general array with many domains the regime of existence may be therefore somewhat 
different. 
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Figure 3: Regime of existence of bound pairs of domain walls for S = 1. The calculation was 
performed with 12800 points and a grid spacing of 0.034. The ramped part is 400 points wide 
on each side. 



3 Application to Two Dimensional Systems: Zig-Zags 

The above results are expected to capture also certain aspects of patterns extended in two 
dimensions. A particularly clear example is given by electro-hydrodynamic convection (EHC) in 
nematic liquid crystals. In the nematic phase the rod- like molecules are predominantly oriented 
in one direction rendering the fluid anisotropic. Depending on parameters, convection arises 
in the form of rolls perpendicular or oblique to that preferred direction. Due to reflection 
symmetry, the oblique rolls can have wave number (Q, P) or (Q, —P). Often domains of oblique 
rolls with opposite orientation are observed to coexist, leading to 'zig-zag' patterns. They are 
characterized by sharp transitions between the domains with the two orientations. 

Quite generally, in isotropic systems, as for instance in Rayleigh-Benard convection (RBC), 
straight roll patterns exhibit a secondary instability which leads to an undulatory deformation 
of the rolls. It arises for small wave numbers and tends to increase the local wave number. 
In the course of the nonlinear evolution the undulations can grow sharper and form zig-zag- 
type patterns which may in fact be stable, as recently shown in simulations of chemical Turing 
patterns p. 

Near onset both systems can be described by suitable Ginzburg-Landau equations which 



differ in tlieir scaling of the spatial variables. In the isotropic case one obtains [23| 

dTA = -{idx + dlfA + \A-\A\^A, (4) 

whereas the anisotropic case yields [^ 



StA = {d\ - iZdxdl + Wdl - d^)A + \A- \A\^A. (5) 

In the anisotropic case the normal/oblique transition occurs aX.W = —ZQ. 

Focussing on solutions which are strictly periodic vn X, A = Ai{Y,T)e^^ , one obtains 
eq.(g)for^i, 

dTAi=D2dlAi-d^Ai + ^A-\Ai\^Ai, (6) 

with 

D2 = Q, T. = \-Q'^ for (4) (7) 

D2 = W + ZQ, S = A-Q2 for(D. (8) 

Note that Y and P in the two-dimensional systems correspond to X and Q in eq.(|l|). Thus, 
the domain structures discussed in sec.|| are one-dimensional analogs of zig-zag patterns. The 
general stability of zig-zag patterns in anisotropic systems has been studied by Pesch and Kramer 
||l7|| and by Bodenschatz et al. [Q] who find islands of stability for the zig-zag states. They are 
presumably related to the discrete set of domains with different widths discussed above. Of 
course, the stability of the zig-zags may be limited by two-dimensional instabilities which may 
preempt some of the transitions discussed here. The interaction of domain walls separating the 
two kinds of oblique rolls has, however, not been discussed in detail. 

The straight (normal) roll state corresponds to a solution with constant Ai. Reducing its 
wave number is equivalent to decreasing D2, which goes through zero at the onset of the zig-zag 
instability. It leads to the growth of undulatory deformations of the rolls which correspond to 
periodic variations of the local wave number in eq.(|^). Based on the simulations of that equation 
one would expect that D2 has to be sufficiently negative for the resulting zig-zag structure to 
persist. Therefore, close to the onset of the zig-zag instability the coarsening dynamics predicted 
by the phase equation dH) will prevail, and the expected final state would be a pattern of straight 
(oblique) rolls of either sign if the phase is not conserved. For deeper quenches into the unstable 
regime]^, however, the undulations are expected to lock into each other due to the oscillatory 
interaction and form an array of domain walls, i.e. an array of zigs and zags, which again could 
be spatially chaotic. 

4 Conclusion 

The competition between patterns differing only in their wave number can lead to complex 
spatial patterns in which the interaction between the walls separating the domains with different 
wave number plays an important roll. Previously, we investigated them for conditions under 



which the total phase, i.e. the number of convection rolls in the system, is conserved [19|. There 
we found a very rich bifurcation structure which originates from interactions with other modes 
arising from the Eckhaus instability. 



^Note |D2/Si ^ 00 at the neutral curve E = A — Q^ 



In experiments the phase is not always conserved. In the thin layers used in electro-convection 
of liquid crystals, for instance, inhomogeneities in the layer thickness are often strong enough that 
convection arises in large patches which are separated by non-convecting regions. In isotropic 
two-dimensional systems focus singularities are expected to arise in the corners of the container 
which also allow rolls to disappear smoothly without strong pinning [^. Here we therefore 
studied domain structures in the presence of spatial ramps. We found that they can exist 
even without phase conservation and we attribute this to the oscillatory behavior of the local 
wave number which should lead to an oscillatory interaction between adjacent domain walls. 
The distance between the walls is discretized accordingly. Close to the onset of the zig-zag 
instability in two-dimensional systems, these oscillations are small; to stabilize zig-zag patterns 
one therefore has to quench deeper into the unstable regime. 

This work has been supported by grants from NSF/AFOSR (DMS-9020289) and DOE (DE- 
FG02-92ER14303). 
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